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Abstract 



We consider the separability of various joint states for A^ qutrits. We derive 

two results: (i) the separability condition for a two-qutrit state that is a 

mixture of the maximally mixed state and a maximally entangled state (such 

^ ' a state is a generalization of the Werner state for two qubits); (ii) upper and 

lower bounds on the size of the neighborhood of separable states surrounding 

f— ^ . the maximally mixed state for N qutrits. 
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*We dedicate this article to Marian Scully and his lifelong interest in the weird and wonderful world 
of quantum mechanics. We understand that the article will appear in the Festschrift celebrating 
Marian's 60th birthday, but surely that 60 is a mistake. Perhaps Marian's passport says he's 60, 
but that's because his passport doesn't know him. He approaches physics with the same curiosity 
and zest he displayed when we met him nearly twenty years ago. Happy 60th, Marian! We draw 
inspiration from your dedication to and enthusiasm for physics. 
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I. INTRODUCTION 

Quantum entanglement provides a powerful physical resource for new kinds of com- 
munication protocols and computation |jl]], which achieve results that cannot be achieved 
classically. Quantum entanglement refers to correlations between the results of measure- 
ments made on component subsystems of a larger physical system that cannot be explained 
in terms of correlations between local classical properties inherent in those same subsys- 
tems. Alternatively, an entangled state cannot be prepared by local operations and local 
measurements on each subsystem. While the nonclassical nature of quantum entanglement 
was recognized by Einstein, Schrodinger, and other pioneers of quantum mechanics, it is only 
recently that a full appreciation of the surprising and complex properties of entanglement 
has begun to emerge. 

We now have a good understanding of entanglement for a pair of qubits 0, a qubit 
being a system with a two-dimensional Hilbert space. Moreover, we have a criterion, the 
partial transposition condition of Peres 0], which determines whether a general state of two 
qubits is entangled and whether a general state of a qubit and a qutrit is entangled @], a 
qutrit being a system whose states live in a three-dimensional Hilbert space. The partial- 
transposition condition fails, however, to provide a criterion for entanglement in other cases, 
where the constituents have higher Hilbert-space dimensions |^,^ or where there are more 
than two constituents. In this paper we consider the entanglement of composite systems 
made of qutrits. In particular, we investigate the separability of various joint states of two 
or more qutrits. A joint state of a composite system is defined to be separable if it can be 
decomposed into a mixture of product states for the constituents. A separable state has no 
quantum entanglement. 

In Sec. 0, we review an efficient operator representation of the states of a qutrit [0, 
particularly the pure states. This operator representation is analogous to the familiar Pauli 
or Bloch-sphere representation of qubit states. The representation is applied to an analysis 
of qutrit entanglement in Sees. |IT| and 0. In Sec. ^, we consider states of two qutrits that 
are a mixture of the maximally mixed state and a maximally entangled state. Such states 
are a generalization of the Werner state for two qubits [^]. We show that the two-qutrit 
mixture is separable if and only if the probability for the maximally entangled state does not 
exceed 1/4. This result is a special case of a general result obtained by Vidal and Tarrach 
H, which gives, for any bipartite system, the separability boundary for any mixture of the 



maximally mixed state with a pure state of the bipartite system. In Sec. |^ we consider 
the separability of A^-qutrit states near the maximally mixed state. We find both lower and 
upper bounds on the size of the neighborhood of separable states around the maximally 
mixed state. Our analysis follows closely that of Braunstein et al. |T^, who analyzed the 
separability of A^-qubit states near the maximally mixed state. 

II. OPERATOR REPRESENTATIONS OF QUTRIT STATES 

In this section we review an efficient operator representation for qutrit states, analo- 
gous to the Pauli or Bloch-sphere representation for qubits. The qutrit representation uses 
the (Hermitian) generators of SU(3) as an operator basis. Our discussion is taken from 



Ref. [0; the reader is referred there for more details (see also Refs. [|TT| and |]T2|)- This op- 
erator representation — and its generalization to higher dimensions — has been used by Acin, 
Latorre, and Pascual [|13| to explore optimal generalized measurements on collections of 
identically prepared spin systems. 

Let |1), |2), and |3) be an orthonormal basis for a qutrit. To denote these states, we 
use Latin letters from the beginning of the alphabet, which take on the values 1, 2, and 3. 
A pure state \'il)) is a superposition of the states \a). Normalization and a choice for the 
arbitrary overall phase allow us to write any pure state as 



\'4)) = e*»sin^cos0|l) + e*^2sin^sin0|2) + cos^|3) . 
One obtains all pure states as the four coordinates vary over the ranges 

< ^,0< 7r/2 , 



f2.11 



(2.2a) 



0<Xi,X2<27r. 



(2.2b) 



To develop an operator representation of qutrit states, we begin with the eight (Hermi- 
tian) generators of SU(3). In the basis |a), these generators have the matrix representations 
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(2.3) 



We use Latin indices from the middle of the alphabet, ranging from 1 to 8, to label these 
generators and related quantities, and we use the summation convention to indicate a sum 
on repeated indices. The generators (|2.3| ) are traceless and satisfy 



AjAfc — - 5jk + djkiXi + ifjki^i ■ 



(2.4) 



The coefficients fjku the structure constants of the Lie group SU(3), are given by the com- 
mutators of the generators and are completely antisymmetric in the three indices. The 
coefficients djki are given by the anti-commutators of the generators and are completely 
symmetric. Values of these coefficients can be found in Ref. [0. 
By supplementing the eight generators with the operator 
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(2.5) 



we obtain a Hermitian operator basis for the space of linear operators in the qutrit Hilbert 
space. This basis is an orthogonal basis, satisfying 

tr(A„A^) = 25^p . (2.6) 

Here and throughout Greek indices run over the values to 8. 
Any density operator can be expanded uniquely as 

P = gC„A„ , (2.7) 

where the (real) expansion coefficients are given by 

Ca = ^tripXa) . (2.8) 



Normalization implies that cq = -\/3/2, so the density operator takes the form 

p = ^(l + cA) = ^(l + c-A). (2.9) 

Here c = CjSj can be regarded as a vector in a real, eight-dimensional vector space, and 
A = Aj-Gj is an operator vector. Using Eq. ( |2.4| ), one finds that 

where the "star" product is defined by 

c-kd = ejdjkiCkdi . (2-11) 

For a pure state, p^ = p, so we must have c ■ c = 3 and c t^t c = \/3c. Defining the 
eight-dimensional unit vector n = c/a/3, we find that any pure state of a qutrit can be 
written as 



\^){ip\ = l{I + V3n-X) = P^, (2.12) 

where n satisfies 



P ,.,.., 3. 



n ■ n = 1 , (2.13a) 

n.^n = n. (2.13b) 

We introduce the notation P^ for a pure state with unit vector n for use later in the paper. 
Equation (p.8|) implies that 



n 



'3 



2 '' " 2 
Applied to the pure state ( |2.1| ), this gives a unit vector 



^tr(pA,) = ^(V'lA,!^) . (2.14) 



n = a/3 I sin^6' sin cos cos(x2 — Xi) ^i + sin^6'sin0cos0sm(x2 — Xi) ^2 



+ - sm^6'(cos^(; 



sin^0) 63 + sin Q cos Q cos cos Xi ©4 



sin d cos d cos sin xi 65 + sin 6* cos d sin cos X2 ee 

1 



sin 6' cos Q sin sin X2 e? + 



2^3 



1 — 3 cos &) eg 



(2.15) 



The pure qutrit states lie on the unit sphere in the eight-dimensional vector space, but 
not all operators on the unit sphere are pure states. For example, of the unit vectors ie^, 
j = 1, . . . , 8, only —eg satisfies the star-product condition (|2.13b| ) and thus is the unit vector 
for a pure state. It is easy to show that the unit vectors that do not satisfy the star-product 
condition do not specify any state, pure or mixed, for they all give operators that have 
negative eigenvalues. The star-product condition ( |2.13b| ) places three constraints on the 
unit vector for a pure state, thus reducing the number of real parameters required to specify 
a pure state from the seven parameters needed to specify an arbitrary eight-dimensional 
unit vector to four parameters, which can be taken to be the four coordinates of Eq. (p.l|) . 

It is useful to notice that 



(^l^')l =tr(pp') 



> 



2n ■ n') . 



(2.16) 



Orthogonal pure states have unit vectors that satisfy n ■ n' = — | and are thus 120° apart. 
The states in an orthonormal basis have unit vectors that lie in a plane at the vertices of an 
equilateral triangle. The density operators that are diagonal in the orthonormal basis are 
the operators on the triangle or in its interior. 

The orthonormal states |a), a = 1, 2, 3, of the original basis have unit vectors n^ that lie 
in the plane spanned by 63 and eg: 



ni 

"2 

ns = -eg 



\/3 1 

—63 + 2 *^8 {Q = T^l'^, = 0, xi and X2 arbitrary). 
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9 = 7r/2, = 7r/2, X\ and X2 arbitrary), 
0, 0, Xi) and X2 arbitrary). 



(2.17a) 

(2.17b) 
(2.17c) 



These unit vectors correspond to the following operator relations for the projectors onto the 
original basis: 



(2.18a) 
(2.18b) 
(2.18c) 
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[l + v^ffA3 + iA,)l 
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l + v^f-^A3 + ^Ag) 


3)(3 = 


4< 


1 - VSAg) . 





For later use, we note here the operator expansions for the transition operators in the original 
basis: 

|1)(2| = i(Ai + 1X2) , |2)(1| = i(Ai - zAs) , (2.18d) 

|1)(3| = i(A4 + ^5) , |3)(1| = i(A4 - tX,) , (2.18e) 

|2)(3| = i(A6 + tXr) , |3)(2| = ^(Ae - tXj) . (2.18f) 

The geodesic distance between pure states, according to the unitarily invariant Fubini- 
Studi metric 0], is given by the Hilbert-space angle between the states, i.e., cossps = 
|('?/'|'?/'')|. For infinitesimally separated pure states, \ip) and \ip') = lip) + \dip), the distance 
gives the Fubini-Studi line element 

dsls = 1 - Mi^')\' = (#l#) - mdi^)\' ■ (2.19) 

We choose to rescale all lengths by a factor of -\/3, giving a line element 

ds"^ = 3rfs|g = dn-dn, (2.20) 

so that the metric on the submanifold of pure states is the metric induced by the natural 
metric on the unit sphere in eight dimensions. In the coordinates of Eq. ( p.l|) , the rescaled 
line element becomes 

ds^ = ?,(de'^ + sin2^rf02 ^ sin^^cos^ll - sin^^cos^) dx\ 

+ sin^^ sinV(l - sin^^ sin^) dxl - 2 sin^^ sin^ cos^0 dxi dx2) ■ (2.21) 

The corresponding unitarily invariant volume element on the space of pure states is 

dVtn = 9 sin'^6' cos 9 sin cos d9 dcj) dxi dx2 , (2.22) 

which gives a total volume 

V = / rffin = (2.23a) 

(notice that the volume relative to the original Fubini-Studi scaling is Vfs = V/9 = 7r^/2). 
We can also show that the components of n satisfy 

dQnnj = 0, (2.23b) 

/ dfinrijUk = -TTT^jk = -3- Sjk ■ (2.23c) 



III. MIXTURES OF MAXIMALLY MIXED 
AND MAXIMALLY ENTANGLED STATES 

In this section we deal with two qutrits, labeled A and B. We consider a class of two- 
qutrit states, specifically mixtures of the maximally mixed state, Mg = ^1 ® 1, with a 
maximally entangled state, which we can choose to be 

1^) = ^(|1) ® |1) + |2) ® |2) + |3) ® |3)) . (3.1) 

Such mixtures have the form 

p, = (l-e)M9 + e|vl>)(vl>|, (3.2) 

where < e < 1. A state of the two qutrits is separable if it can be written as an ensemble 
of product states. In this section we show that the state ( |3.2| ) is separable if and only if 
e < 1/4. 

In analogy to Eq. ( p. 71) , any state p of two qutrits can be expanded uniquely as 

P = -CapK ® A^ , (3.3) 

y 

where the expansion coefficients are given by 
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Cap = -rtripXa ® A;3) • (3.4) 



Normalization requires that cqo = 3/2. Using the relations ( p.l8|) , we can find the operator 
expansion for the maximally entangled state ([3.11): 



|v[/)(vl/| = ^|a)(6|®|a)(6| 

a,b 

1 / 3 

= - 1 O 1 + -(Ai O Ai - A2 O A2 + A3 ® A3 

+ A4 ® A4 - A5 ® As + Ae ® Ae - A7 ® A7 + As ® As) I . (3.5) 

Hence the operator expansion for the mixed state ( p.2|) is 
1 / "^ 

P, = - I 1 ® 1 + y (Ai ® Al - A2 ® A2 + A3 ® A3 

+ A4 ® A4 - As ® A5 + As ® As - A7 ® A7 + As ® As) J , (3.6) 

from which we can read off the expansion coefficients ( p.4|) for the state p^: 

coj = Cjo = , Cjk = for j ^ k, (3.7a) 

Q 

cii = -C22 = C33 = C44 = -css = C66 = -C77 = css = y • (3.7b) 
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The product pure states for two qutrits, P^^ ® Png , constitute an overcomplete operator 
basis. Thus we can expand any two-qutrit density operator in terms of them: 



P 



I dVL^^ dn^B w{nA, nB)PnA ® Pub ■ (3-8) 



Because of the overcompleteness, the expansion function w{jia^'^b) is not unique. Notice 
that the expansion coefficients Cafs of Eq. (|3.4|) can be written as integrals over the expansion 
function, 

CafS = S dO^nA d^ng willA, n^) (nA)a(nB)/3 , (3.9) 

where uq = l/-\/2 and hj = rij. 

A two-qutrit state p is separable if there exists an expansion function that is everywhere 
nonnegative. For a separable qutrit state, the expansion function w(n^, n^) can be thought 
of as a normalized classical probability distribution for the unit vectors n^ and n^, and 
the integral for Cq/j in Eq. ( p.9| ) can be interpreted as a classical expectation value over this 
distribution, i.e., 

Ca/3 = 3E[{nA)a{nB)p] . (3.10) 

If the state p^ is separable, we have from Eqs. ( p.7b| ) and (|3.10|) that 



e _ 1 
2 ~ s'""^'^'' 



E[{nA),{nB),]\ < ^{E[inAr,]+E[{nB)',]) . (3.11) 



Adding over the eight values of j gives 

4e < ^{E[nA ■ n^] + ^[n^ ■ n^]) = 1 . (3.12) 

We can conclude that if pe is separable, then e < 1/4. 

To prove the converse, we need to construct a product ensemble for e = 1/4. For each 
pair a, b, with b > a, define four pure product states 

where z takes on the values ±1 and ±i, so that = J2z ^ = J2z ^"^ and 4 = J2z kP- It is easy 
to show that an ensemble of the twelve states l^^"'*-*), all states contributing with the same 
probability, gives the state Pe=i/4, i.e., 

^ E E |$^''))($^'')| = Im, + \\^){^\ , (3.14) 

a,b ^ 
b>a 

thus concluding the proof that p^ is separable if and only if e = 1/4. 

This result should be contrasted with that for the corresponding qubit state, a Werner 
state ip, which is separable if and only if e < 1/3 [|T5|. This indicates that maximally 
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entangled states of two qutrits are more entangled than maximally entangled states of two 
qubits. 

The separability boundary obtained in this section is a special case of a general result 
obtained by Vidal and Tarrach |^, who showed the following: for a pair of systems, one 
with dimension di and the other with dimension d2, a mixture, (1 — e)Md^d2 + ^l^)(^|) 
of the maximally mixed state Mci^d2 with any pure state \^) is separable if and only if 
e < {1 + did2aia2)~^ 1 where a\ and a\ are the two largest eigenvalues of the marginal density 
operators obtained from |\E')(\1'|. The reason for presenting the much more limited result of 
this section is, first, that the proof of necessity, culminating in Eq. ( ^.121 ), has a nice physical 
interpretation in terms of the correlation coefficients E[{nA)j{nB)j\ and, second, that the 
ensemble (|3.14|) is different from the one used by Vidal and Tarrach. 



IV. SEPARABILITY OF STATES NEAR THE MAXIMALLY MIXED STATE 

This section deals with A^-qutrit states of the form 

p, = (1 - t)M^N + epi , (4.1) 

where M^n = I®- ■ -(8)1/3^ is the maximally mixed state for N qutrits and pi is any A^-qutrit 
density operator. In this section we establish upper and lower bounds on the size of the 
neighborhood of separable states surrounding the maximally mixed state. In particular, we 
show, first, that for e < (1 + 3^^~^)~^, all states of the form ( ^.1| ) are separable and, second, 
that for e > (1 + 3^/^)~^, there are states of the form ( [4.1|) that are not separable. The 
approach we take in this section follows slavishly the corresponding qubit analysis presented 
by Braunstein et al. |[T0| . 

Since the product pure states form an overcomplete operator basis, any A^-qutrit state 
can be expanded as 



P 



j dn^, ■ ■ ■ dn^^ w(ni, . . . , njv)Pni ® ■ ■ ■ ® i^n^ • (4.2) 



The expansion function w{ni, . . . , n^v) is not unique, because of the overcompleteness of the 
pure product states. An A^-qutrit state is separable if there exists an expansion function 
that is everywhere nonnegative. 

We can find a particular expansion function in the following way. Any A^-qutrit density 
operator can be expanded uniquely as 

P= ^Car...aj,K,®---®KM , (4-3) 

with the expansion coefficients given by 



'--ai...aiv 



/3\ 

= i^-j tr(pA„, ®---®A,^). (4.4) 

Using Eqs. (|2.23| ), we can write 

16 /• _ 

K = nr ^ J dQn riaPn , (4.5) 



where the barred components are defined by uq = l/4-\/2 and rij = rij. Substituting this 
expression for Aq, into Eq. (f4.3|), we find that one choice for the expansion function is 



/ 16 

N 



N 



Wp[ni,...,nN) = I Q /o^2 I Cai...a^(nijc.i---ln7vja^ 



2 

' tr(p(l + 4V3ni-A)®---®(l + 4v^njvA)) . (4.6) 



.97r2 

Now consider the product operator in Eq. ( |4.6|) . The pure-state density operator P^ = 
|(1 + v^n ■ A) has eigenvalues 1, 0, and 0. This imphes that n ■ A has eigenvalues 2/\/3, 
— 1/-\/3, and — 1/-\/3 and that each term in the product operator, 1+4-^/3 n- A, has eigenvalues 
9, -3, and -3. Hence the smallest eigenvalue of the product operator is 9^^^ (—3) = — 3^^"^. 
The result is a lower bound on the expansion function (|4.6|): 



/smallest eigenvalue of\ / 2 \ ^2N-i 



^^ ' ' ^ - V97rV V product operator / V97r2/ ^ ^ 

We can use this lower bound to place a similar lower bound on the expansion function 
for a state of the form ( [4.1[ ). Since the expansion function for the maximally mixed state 
M^N is the uniform distribution {2/9tt'^)^ , we have 

w,Xnu...,n^) = {l-e)[^) +^^P. >(^) (1-6(1 + 3^^-1)). (4.8) 



We conclude that if 



1 

1 + 32 



^ < TT^^^Fn ' (4-9) 



Wp^ is nonnegative, and thus the qutrit state pe of Eq. (f4.1|) is separable. This establishes 
a lower bound on the size of the separable neighborhood surrounding the maximally mixed 
state. 

We turn now to obtaining an upper bound on the size of the separable neighborhood. 
Consider two particles, each with spin (3^^^ — l)/2 (A^ even) and thus each having a {d = 
3^/2)-dimensional Hilbert space. We can think of each particle as being an aggregate of 
N/2 spin-1 particles (qutrits). We consider the following joint density operator for the two 
particles (or of the A^ qutrits), 

p, = (l-e)M,2+e|</))(</)|, (4.10) 

where M^^ = 1/(P' is the maximally mixed state in (P = 3^ dimensions, and 

10) = — (|1) ® |1) + |2) ® |2) + • ■ ■ + \d) ® \d)) (4.11) 

is a maximally entangled state for the two particles. 

Now project each particle onto the subspace spanned by |1), |2), and |3). For each 
particle this subspace can be thought of as the Hilbert space of a single qutrit; thus the 
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projection leaves a joint state of two qutrits. The projection is effected by the projection 
operator 

3 3 

n=^|a)(a|®^|6)(6|, (4.12) 

a=l b=l 

which is the unit operator in the projected two-qutrit space, i.e., 11 = 1 (8> 1 = 9Mg. The 
normahzed state after projection is 

where 

A = tr(p,n) = (1 - ^)^ + f = ^(l + <d/^ - 1)) (4-14) 

is a normahzation factor, |\[^) is the maximally entangled state of two qutrits given in 
Eq. Q, and 

A 1 + e(d/3 - 1) ^ ^ 



The projected state p^ is the state ( p.2|) considered in Sec. |T|, a mixture of the maximally 



mixed state for two qutrits, Mg, and the maximally entangled state |\1/). As shown in Sec. |T|, 
this state is nonseparable for e' > 1/4, which is equivalent to 

Moreover, since the local projections on the two particles cannot create entanglement from 
a separable state, we can conclude that the state ( |4.1CI| ) of A^ qutrits is nonseparable under 
the same conditions. This result establishes an upper bound, scaling like 3~^/^, on the size 
of the separable neighborhood around the maximally mixed state. 
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